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Wc study reconnections of quantum vortices by numerically solving the governing Gross-Pitaevskii equation. 
I We find that the minimum distance between vortices scales differently with time before and after the vor- 

04 ' tex reconnection. We also compute vortex reconnections using the Biot-Savart law for vortex filaments of 

^ ' infinitesimal thickness, and find that, in this model, reconnection are time-symmetric. We argue that the 

Q . likely cause of the difference between the Gross-Pitaevskii model and the Biot-Savart model is the intense 

' rarefaction wave which is radiated away from a Gross-Pitaeveskii reconnection. Finally we compare our re- 

sults to experimental observations in superfluid helium, and discuss the different length scales probed by the 
' two models and by experiments. 

PACS numbers: 
^ I 47.32.C- (vortex dynamics) 

fH , 67. 30. he (vortices in superfluid helium) 

03.75.Lm (vortices in Bose Einstein condensates) 



I. INTRODUCTION 



The importance of vortex reconnections in turbulence^ cannot be understated. Reconnections randomize the velocity 
field, play a role in the energy cascade, contribute to the fine-scale mixing, enhancing diffusion^Tii, and are the dominant 
^ ' mechanism of jet noise generation. If the axes of tubular vortex structures are interpreted as the skeleton of turbulence, 
then the knottcdness of the axes characterizes the turbulence's topology, and vortex reconnections are the critical 
events which change this topolog}*^. This idealized picture becomes physical reality if one moves from ordinary viscous 
^ ' fluids to quantum fluids^ such as superfluid liquid helium ('^He-B and ^Hc) and atomic Bose-Einstein condensates. 

, In these superfluid systems, quantum mechanics constrains any rotational motion to vortex lines around which the 
^\ ' circulation is fixed by the condition 

(N 

I u dr= - = K, (1) 
O ■ /c 
<N . 

where u is the velocity, C is a closed integration path around the vortex axis, h is Planck's constant, m is the mass 
of the relevant boson (a helium atom in the case of ''He, a Cooper pair in the case of '^He-B), and k is the quantum 
, , , of circulation. Another constraint of quantum mechanics is the small vortex core, which has fixed radius (10~^ cm 
rN ' in '*He, 10~^ cm in ^He-B) and is orders of magnitude smaller than the average distance between vortices in typical 
I experiments; because of these constraints, no intensification or diminution of vorticity through stretching of the vortex 
" " " core is possible in quantum fluids. 

Vortex reconnections of individual quantum vortex lines are discrete, dramatic events, which have been recently 
visualized in the laboratory^. They are the key to understanding quantum turbulence^, a disordered state of vortex 
lines which is easily created by stirring liquid helium and atomic condensates. Compared to quantum reconnections, 
viscous reconnections in ordinary (classical) fluids are not complete events: classical vorticity is continuous, not 
discrete, and parts of the initial vortical tubes can be left behind as vortex threads, which then undergo successive 
reconnections (the cascade and mixing scenarios), as newly formed vortex bridges recoil from each other by self- 
advectionS. Another important difference is that classical reconnections are dissipative events: viscous forces turn 
part of the fluid's energy into heat, whereas in a superfluid the viscosity is zero, and the fluid's total energy is 
conserved. 

It can be argued that, because of the utter simplicity of quantum vortices (zero viscosity, fixed circulation, small 
core size), quantum vortex reconnections are not only important phenomena of low temperature physics and of atomic 
physics, but are also relevant to our general understanding of fluid phenomena as toy models of Euler dynamics. 
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Although a rather large number of studies addressed the viscous rcconncction problem, both the underlying mecha- 
nism and various scaling relationships have remained elusive^. The first analytical work on classical vortex reconnec- 
tions goes back to Crowi^, who studied the instability of a pair of counter-rotating vortex tubes shed from the wing 
tips of an airplane. However, the systematic study of vortex reconnections began with observations and laboratory 
experiments of the simplest and most fundamental interaction of two colliding vortex ringsiir— . With the rapid 
development of supercomputers, direct numerical study of vortex reconnection became possible. Ashurst & MeironM 
numerically solved the incompressible Navier-Stokes equations in the region of closest approach of two vortex rings 
by providing an initial condition generated using the Biot-Savart (BS) modeli^ of vortex filaments. At the same time, 
Pumir & Kerr— performed numerical simulations of interacting vortex tubes. Many other studies followed (for a 
review of the extensive work up to 1994 see the review of Kida & Takaoka^), but the literature on classical vortex 
reconnections seemed to fade after 1994i2r— . The recent work by Hussain & Duraisamy^ renewed the interest in the 
mechanics of viscous vortex reconnection. Their study focused on the direct, high resolution numerical simulation 
of the incompressible Navier-Stokes equations over a wide range of vortex Reynolds numbers. Re, for two perturbed 
anti-parallel vortex tubes. They found that the minimum distance S between the vortex tube centroids scales as 
S{t) ^ {to — t)^/^ before the reconnection and as S{t) ^ {t ~ to)^ after the reconnection, where t is time and to is the 
instant of smallest separation between the vortex centroids. 

The literature concerned with quantum reconnections is more limited. The possibility of reconnections was 
first raised by Fcynman^S. in his pioneering work on the applications of quantum mechanics to superfluid helium. 
Schwara^^^ realized that vortex reconnections are necessary to account for quantum turbulence. He modelled quan- 
tum vortices as classical vortex filaments and proposed the Local Induction Approximation (LIA^^) as a practical 
alternative to the exact (but CPU-intensive) Biot-Savart law for the numerical study of superfluid vortex dynamics. 

A few years later, Koplik & Levine^^ performed the first numerical simulation of quantum vortex reconnections by 
solving the Gross-Pitaevskii equation (GPE) which governs the motion of a Bose-Einstein condensate and is used as a 
model of superfluid helium. They found that if two vortices are nearly anti-parallel when the large-scale fluid motion 
brings them together, they reconnect, thus confirming the conjecture of Feynman and Schwarz. 

Soon after, using the vortex filament model of Schwarz, de Waele & Aarts^i numerically integrated the Biot-Savart 
equation for the ideal-fluid velocity field and claimed the existence of a universal route to reconnection for all kind 
of initial vortex-antivortex arrangements: their calculations showed that, when vortices approach each other, they 
always form the same pyramidal cusp. They measured the minimum distance between vortices, 6, during the approach 
to reconnection as a function of time t, and found that 



S{t)^^-{to~t), (2) 

where to is the time of reconnection. Equation ([2]) is consistent with the dimensional argument that, if the only 
relevant parameter in reconnection dynamics is the quantum of circulation, then 

5it)^Ai^\to~t\r, (3) 

where a = 0.5 and A is a dimensionless constant of order unity. 

Leadbeater et al^ used the GPE model to study reconnections of vortex rings launched against each other, and 
discovered that a sound wave (in the form of a short rarefaction pulse) is emitted at a reconnection event. The wave 
turns part of the initial kinetic energy of the vortices into acoustic energy which is radiated away. Acoustic energy 
is also created by vortex acceleration^Si^S. The effects are clearly important to make sense of the observed decay of 
quantum turbulence at very low temperatures^. 

Individual quantum vortex reconnections were first observed by Paoletti et al^ by analyzing the trajectories of 
solid hydrogen tracers in superfluid ^He. They verified that the scaling with a = 0.5 holds before and after the 
reconnections, that is to say quantum reconnections are statistically time reversible, unlike classical reconnections^. 
They concluded that the expression 

Sit)^AiK\to-t\y^'{l + c\to~t\), (4) 

is the best fit to their data, with wide distribution of values centred around A sa 1.1 (larger than A = (27r)~^/^ « 0.4 
found by de Waele & Aarta^i) and c sa 0. 

The universality of the route to reconnections was questioned by Tebbs et al^ who performed a series of numerical 
simulations of quantum reconnections using the GPE. They reproduced the pyramidal shape of the vortex lines 
observed by de Waele & Aarts^ for the initial configuration used by these authors, but did not observe the same 
shape for other configurations. They also confirmed the scaling (jS]), again in agreement with de Waele & Aarts— , 
but did not measure the time evolution of S after the reconnection. 
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Kursa et alM- employed LIA, Biot-Savart and GPE simulations, and showed that a single reconnection of two 
almost anti-parallel quantum vortices can lead to the creation of a cascade of vortex rings, provided that the angle 
between the vortices is sufficiently small. 

Kerr—, by means of numerical simulation of the GPE, investigated the reconnection between a pair of perturbed 
anti-parallel quantum vortices. He argued that kinetic energy is converted into interaction energy and eventually local 
kinetic energy depletion that is similar to energy decay in a classical fluid, even though the governing equations are 
Hamiltonian and energy conserving. 

The present study aims at characterizing vortex reconnection in quantum fluids by performing direct numerical 
simulations of the three-dimensional GPE in different vortex configurations. Our goal is to extract the minimum 
distance S between vortices as a function of time, both before and after the reconnection, compare results against 
classical reconnections, quantum reconnections computed with the Biot-Savart law, and experiments in superfluid 
helium. 

The plan of the paper is the following. In section|lI]we introduce the governing GPE, cast it in dimensionless form, 
and present the straight vortex solution. In section Hill we describe the initial condition of our numerical calculations 
in terms of the initial angle /3 between the vortex lines, and present computed vortex reconnections, paying attention 
to the minimum distance S between vortices before and after the reconnection, and the sound wave which is generated. 
We also perform reconnections using the Biot-Savart model. In section ITVl we compare GPE reconnections and Biot- 
Savart reconnections with previous work and experiments, and draw the conclusions. The numerical method which 
we use to solve the GPE is described in the Appendix. 



II. MODEL 

The governing equation for a weakly-interacting Bose-Einstein condensate is the GPE)2^ 



ift^ =-^VV + T^o|V'l'V'-£^o^, (5) 

where ■(/'(x,t) the macroscopic single-particle wave function for N bosons of mass m at position x and time t, Vq is 
the strength of the repulsive interaction between the bosons, Eq is the chemical potential (the energy increase upon 
adding a boson), h = /i/(27r), and the normalization condition is 



J |V'pdx = TV. (6) 

To study atomic condensates rather than liquid helium, a term of the form Vtrap^, where Vtrap is a suitable trapping 
potential (usually harmonic), is added to the right hand side of equation [5l Hereafter we shall not consider such term, 
but limit our work to homogeneous condensates. When applying the GPE to superfluid helium we must remember 
that helium is a liquid, not a weakly interacting gas, so the GPE model is more qualitative than quantitative. For 
example, the dispersion relation of small perturbations from the uniform solution ^JJao = \/Eo/Vo of equation [5] is 



2 _ ^^^'"^ , ^0,2 



k', (7) 



where ui is the angular velocity, k is the wavenumber and c — \J EqItu is the speed of sound. Note that for fc << I we 
have Lo K, ck (phonons), and for k >> 1 we have ui « hk"^ /(2m) (free particles), without the roton minimum which is 
characteristic of superfluid helium^. Another shortcoming of the GPE is the vortex core: more sophisticated models of 
the helium vortex core cxist^^, but are not practical for the study of complex dynamics such as vortex reconnections. 
By applying the transformation 



X —= X, t ——t, 

^/2^ 2Eo' V 1^0 



where tpoa = \/EqIVq is the uniform solution at rest in an infinite domain, we cast the GPE in the following 
dimensionless form 
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^ = Avv+^(i-i^n^. (9) 

The quantity Co = fi/^/2mEo is called the healing length or coherenee length. It is the typical length scale over which 
the wave function bends, and therefore determines the vortex core radius (see next section) and the thickness of any 
superfluid boundary layer in the presence of a wall. 

The numerical method to numericaly solve the GPE is described in the Appendix. Figure [T] (left) shows that during 
the time evolution the Hamiltonian energy, defined by 



(10) 



is sufficiently well conserved: the numerical error is less than one part in 10^. The fluid dynamics interpretation of 
the GPE arises from the Madelung transformation 



ij = y^c'', (11) 



which yields the following equations 



dp _^ djpuj) 

dt dxi 



0, (12) 



dui dui \ dp dr, 



M^+^.ie =-7e + ld^' 1-2,3 (13) 



dt duj J dxi dx 



(written in tensorial notation), where density and velocity are 

p=|^|2, u = V5, (14) 

and 

^ = T "- = 4^a^' ('^^ 

are the pressure and the so-called quantum stress. It is easy to verify that the quantum stress term at the right hand 
side of in equation ((T^ is negligible compared to the pressure term at length scales larger than unity (the coherence 
length). Since, as we shall see, the vortex core radius is of the order of the coherence length, we conclude that, at scales 
larger than the vortex core, the GPE, expressed by equations and reduces to the classical (compressible) 
Euler equations. Incompressible Euler dynamics is achieved in the further limit of small velocity (compared to the 
sound speed c = y^Eo/m) at constant density (again, away from the vortex cores). 

We seek a two-dimensional solution of equation ^ that represents a straight vortex centred at the origin. It is 
well-known that the classical two-dimensional Euler vortex of circulation P has azimuthal velocity ug ~ r/(27rr) where 
r = \fx^\-y^ is the radius and = arctan(j//a;) is the azimuthal angle. The Cartesian components of the velocity 
are thus Ux = — sin0 = —Tyj (27rr^) and Uy ~ cos d = Txjr^ . Therefore u = (u^^Uy) = (P/(27r)) Vf?. This shows 
that the velocity field is solenoidal (V • u = 0), that the quantum mechanical phase, 5, is simply the azimuthal angle 
0, and that the quantum of circulation, in our dimensionless units, is equal to 27r. In steady conditions, the continuity 
equation ensures that V • (pu) = 0, hence u • Vp = 0, which means that Vp • = 0. The solution p = const has 
infinite energy and must be rejected. The other possibility is that Vp _L V0. Since the cylindrical coordinates r 
and d are perpendicular to each other, V0 is parallel to r, hence Vp is perpendicular to r. Therefore, for a steady 
two-dimensional and divergence- free velocity field, we must have p = p(r). Seeking a solution which represents a 
vortex centred at the origin, we set -0 ~ ^/pc^^ = /('')g'^ where /(r) is a function to be determined. By imposing that 
■0 is the steady solution of equation we find that /(r) satisfies the equation 
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r + -+f(i-f-^)=o (16) 

with boundary conditions /(O) — 0, /(oo) ~ 1. The equation can be integrated numerically, but, to make the 
computation of the initial condition faster, we look for a Pade approximation R[r) to /(r) of the form 

p. ^ Er=o V}-^' _ Po + Pir + P2r^ + ■ ■ ■ + Prnr" 

1 + 2^/c=i ^fc'' 1 + 91'' + 92?"^ H 1- gn'^" 

which agrees with /(r) at the origin to the highest possible order, i.e. /(O) i?(0), /'(O) = i?'(0), /"(O) = i?"(0), 
. . . , /('"+") (0) = i?'^™+"'(0). In our case p > 0, and the Pade approximation can be limited to^ 



- :yxf\ . (18) 

with 

11 5-32ci / lA 

^^^32' ^^ = 48- 192ci ' ^^-^^(.^^-ij- 

The quantum vortex core is thus a hole of radius of the order of ao around which the quantum mechanical phase 
changes by 27r. Figure [T] (right) compares the radial profile of the density near the vortex axis computed numerically 
and with the Pade approximation. 



III. RESULTS 



A. GPE initial condition 



Our initial condition consists of two straight vortex lines Vi and V2 as in figure [2] (the arrow denotes the direction 
of the vorticity), which intersect the y-axis respectively at the points C\ and Ci- We call tti the y-z plane, ni the 
plane perpendicular to tti and passing through C2, and s the intersection between tti and t^i. We call /? the angle 
formed by the directions v'j^ and V2, where v'j^ is the projection of vi onto ■712- In the case of two parallel vortices we 
have /3 = 0; in the case of two anti-parallel vortices we have /3 = tt. Perpendicular intersections occur for /3 = 7r/2 
and /3 = —it/ 2 (the latter is show at the right of figure [5]). 



B. GPE reconnection of anti-parallel vortices 



Our first numerical experiment is concerned with the reconnection of anti-parallel vortex lines (reconnecting angle 
/3 = tt). The size of the computational box is —30 < x,y, z < 30. At time t ~ the vortex lines are located at position 
(xQiyo) — (10; ±3). In order to make sure that the reconnection occurs in the centre of the computational domain, 
we impose a slight initial perturbation to the vortex lines, of the form A [cos(27r(z,nax — Zj^in) / ^)]^ ■ Figure ([3]) shows 
some snapshots of the time-evolution. It is apparent that the two vortices move as a pair in the x-direction. The 
slight initial curvature enhances the Crow instability and the vortices approach each other. The reconnection results 
in the formation of two U-shaped vortex lines, which bend and move apart from each other. 

Figure ([¥]) shows the minimum distance 6 between reconnecting vortices as a function of time t. The algorithm 
to compute 6 searches for vortex lines, starting from the boundaries, following the minimum of the density p. Once 
the vortex lines are retrieved, they are put in parametric form as curves in E'^, i.e. r(^) = {x{^),y{S,), z{^)), and S 
is computed. Our main finding is that 6{t) behaves differently before and after the reconnection. In this respect, 
quantum vortex reconnections are therefore similar to viscous reconnections. Indeed figure (jS]) resembles very closely 
figure (3b) of Hussain & Duraisamy^. 
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C. GPE reconnection of vortices at difTerent initial angles 

We perform several other numerical calculations of vortex reconnections to check whether the time dependence 
of S{t) depends on the angle /? of the initial condition. Figure ([S]) shows the time evolution of two initially straight 
vortices set at the angle f3 = 37r/4. Again, we observe that the vortices first approach each other, then move away after 
the reconnection. A very similar behavior characterize vortices initially set at /3 = 7r/2 (orthogonal reconnection), 
whose time-evolution is shown in figure 

Figure ([7]) summarizes the results obtained for different values of /3. The minimum distance between the vortices, S, 
is reported separately before and after reconnection as a function of |t — to|j where to is the time at which the vortex 
reconnection takes place. Power laws (black solid lines) of the form 

6^A\t~tor, (20) 

(where d and t — Iq are dimensionlcss) arc super-imposed to fit the numerical data; the fitting coefficients A and a 
are reported in table HI From Figure ^ and table (jl]) we conclude that GPE reconnections are not time-symmetric: 
the average values of a is a = 0.39 before the reconnection and a = 0.68 after the reconnection. If we average 
all values, before and after, we obtain a = 0.53, which is in fair agreement with a = 0.5 predicted by the scaling 
argument ([3]). Table also shows that the average values of A are A = 1.29 before the reconnection and A = 1.54 
after the reconnection. If we set a = 0.5 and return to dimensional variables, we obtain A = 1.29/A/27r = 0.52 and 
1.54/-\/27r = 0.62, which are about half of the value A ss 1.1 found by Paoletti et al^ in their experiment. 

As a check, we repeat all numerical simulations for opposite angles (e.g. /3 = (3 = ±|7r, etc.) obtaining the 
same temporal dependence of S upon time t. 

D. GPE reconnection wave 

Our calculation confirms the finding of Leadbeater et al^, that a vortex reconnection creates a wave. First we 
consider the reconnection of anti-parallel vortices (/? = 7r/2). By extracting the iso-surfaccs at quite large level, as 
done in figure (|8]) for p = 0.94, one notices the formation of a mushroom-shaped pressure (density) rarefaction wave 
generated by the reconnection. The wave becomes shallower as it moves away from the vortices. The bottom plot 
of figure ([8]) shows the vortex lines together with the pressure wave obtained by replacing the vortex tubes from the 
isosurface at p = 0.94 with the isosurface at p = 0.2. The footprint of the wave is particularly visible in contours of p 
on the plane y = 0, as shown in figure (|9]). 

The pressure wave is clearly visible for relatively small angles between vortices, /? < 7r/2, whereas for larger angles 
it becomes difficult to clearly track it and visualize it. Figure (fTO]) shows the time evolution of the mushroom-shaped 
pressure wave ejected after reconnection for (5 ~ 7tt/8 (isosurfaccs at p = 0.94). 

E. GPE reconnections and vortex rings 

As we have mentioned in SectionU secondary generation of vortex rings following a reconnection event was observed 
in the numerical simulations of Kursa et al.— and Kerr—. Kursa et al. also studied how the emission of vortex rings 
depends on the initial angle between the vortices (almost antiparallel configurations favour the generation of vortex 
rings following the Crow instability). We do not investigate further the generation of vortex rings, since it was already 
studied in detail in cited works. We only remark that if we make our computational box longer in the z direction, 
vortex rings generation becomes visible, as shown in figure[TTJ As for the physical significance of this effect, we notice 
that, according to a recent study of Baggaley et al4i, the distribution of reconnecting angles /3 depends on the nature 
of the quantum turbulence. Quantum turbulence generated with grids or propellers seems classical in nature (for 
example the kinetic energy is distributed on the length scales according to the classical Kolmogorov fc^^/"^ law where 
k is the wavenumber), and contains coherent bundles of vortices which induce reconnections at small angles /?. On 
the contrary, quantum turbulence generated thermally (e.g. counterflow turbulence) is spatially more random, and 
reconnections tend to be antiparallel (/? « tt). 

F. Biot-Savart reconnections 

It is instructive to compare reconnections computed with the GPE with reconnections computed with the Biot- 
Savart law. The latter, which is widely used to study quantum turbulence, approximates vortex lines as space curves 
s = s(^,t) of infinitesimal thickness which move according to 
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ds 

Utt 



i 




X dr, 



(21) 



where ^ is arc length and the hne integral extends over the entire vortex configuration. Equation (j2ip expresses 
incompressible Euler dynamics in integral forroi^. Physically, it assumes that the density of the fluid is constant (zero 
Mach number limit) and that the vortex core is much smaller than any other length scale in the flow (a small parameter 
must be introduced to de-singularise the integral). Vortex reconnections are forbidden by Euler dynamics, therefore, 
when applying equation (|2ip to superfluid helium, we must supplement it with an algorithmic reconnection procedure 
which changes the topology of two vortex filaments when the distance between them is less than a prescribed cutoff 
value, as first explained by Schwar a^^i^^ . The numerical techniques which we use to compute Biot-Savart evolution are 
described in our previous papers^i^. Here it suffices to say that the vortex filaments are discretized into a variable 
number of points Sj (j = 1, • • ■ N), holding their relative distance approximately between and A^/2 where 
represents the prescribed numerical resolution. The reconnection algorithm which we use, which is triggered when the 
vortex separation is closer than A^/2, has been already described in detail^ and compared to other algorithms used 
in the liquid helium literature. It must be stressed that, unlike some of our recent work on quantum turbulence^, the 
results which we present here do not use a tree-algorithm^ to approximate and speed up the calculation of Biot-Savart 
integrals. 

We perform our calculations in an open domain with numerical resolution A^ = 0.0005 cm. The initial condition 
of the first numerical calculation which we present consists of two vortex rings of the same polarity and radius 
R = 0.0477 cm set parallel to each other, side-by-side on the ccz-plane and travelling in the y-direction, initially at 
distance Ax = 0.002 cm from each other. The initial number of discretization points for the two rings is TV = 1600. 
The same initial condition was used by de Waele & Aarts^l. Figure shows the time evolution. Note that the 
resulting vortex reconnection is locally anti-parallel (/3 = tt) . The initial condition of the second numerical calculation 
consists of the same two rings, but initially set perpendicular to each other. The time evolution is shown in figure (jl3p . 
Note that the resulting reconnection is locally orthogonal (/? = ti'/S). 

The minimum distance between vortices, 6, for both parallel and perpendicular rings, is shown in figure p^ . It is 
apparent that the temporal scaling is time symmetric before and after the reconnection, in agreement with equation ^ 
with a = 0.5 and with the results of de Waele & Aarts^i. The time-symmetry of Biot-Savart reconnections contrast the 
time-asymmetry of GPE reconnections showed in the previous sections and of classical reconnections^. We tested the 
dependence of this result on the reconnection algorithm used^, and found no difference in the scaling. This is perhaps 
not surprising as a change of the reconnection algorithm would imply a change of S of the order of A^ « 10^** cm 
only, whilst we measure the evolution of 5 up to distances of the order 10"^ cm. 

The coefficient A however is not the same in all cases (although, for the approach of parallel rings, it is in fair 
agreement with de Waele & Aarts^l). A similar spread was observed by Tsubota & Adachi^. The speed at which 
the vortex lines move away from each other after the reconnection is faster than the speed at which they approach 
each other; this effect is also visible in figures (HH), and qualitatively consistent with the findings obtained with the 
GPE, see figure 0. 

IV. CONCLUSION 

Hussain & Duraisamy^ have shown that, in ordinary incompressible viscous fluids, the minimum separation 6 
between reconnecting vortex tubes behaves differently before {S{t) ~ {to — t)^^"^) and after {6{t) ~ (< — ^o)^) the 
reconnection at t = t^. By solving the GPE we find a similar time asymmetry, although with different power laws: 
5{t) ~ (to — t)'^'^ and 5{t) ~ (t — ^o)"'^ respectively, independently of the initial angle between the vortex lines. On 
the contrary, by solving the Biot-Savart equation, we find that the scaling is time symmetric, with S{t) ~ \to — t\^^'^ 
for both t < to and t > to- 

What causes the difference? The main difference between the GPE model and the Biot-Savart model is that the 
former is compressible and the latter is not. Clearly the rarefaction wave which is generated at the GPE reconnection 
breaks the time symmetry, transforming^^ some of the kinetic energy of the vortices into acoustic energy which is 
radiated to infinity (in analogy with the viscous dissipation of kinetic energy at Navier-Stokes reconnections). The 
fact that GPE reconnections do not follow the power law ^ predicted by the simple dimensional argument is not 
surprising: at the very small scales explored by solving the GPE, other parameters besides the quantum of circulation 
may be relevant in determing i5, for example the coherence length. 

It is also clear that the Biot-Savart model and the GPE model probe different length scales. As remarked in 
Section ini the GPE should converge to incompressible Euler behaviour in the limit v/c << 1 of small velocity v 
compared to the speed of sound c. The velocity of a vortex strand of local curvature R is approximately v « k/J?, 
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hence we expect to recover the Biot-Savart law for R >> k/c. However, as noticed earher, the GPE is only a 
qualitative model of helium. In estimating k/c we should not use the observed value c = 238 m/s, but rather the 
value of c which arises from the GPE itself, consistently with the vortex core size resulting from the GPE. The first 
step is to identify the coherence length Co- Experiments with ions and vortex rings by Rayfield and Reif suggest^^i^S 
that the radius of the vortex core is og = 1.3 x 10~* cm. Figure [T] (right) shows that the density raises from zero to 
half of its value at infinity at distance 1.5^0 (1-5 dimensionless units) from the axis of the vortex. Taking this distance 
as the (arbitrary but reasonable) value of the vortex core radius oq, wc have Co — 0.87 x 10~* cm. The half-size of 
our computational box (30 coherence length units) thus corresponds to 26 x 10"^ cm. The second step is to find the 
sound speed in the GPE evolution. From c = yjEofm and Co = hj^/lmEi^, we have c = h/ {■\/2'mQQ) = 129 m/s, 
where m = 6.64 x 10~^^ g is the mass of one helium atom. Wc conclude that GPE evolution should become similar 
to Biot-Savart evolution for R >> k/c ^ 8 x 10~® cm, that is to say for i? >> 9Co (radius of curvature much greater 
than 9 dimensionless units). It is apparent in figures [31 [S] and [B] that the radius of curvature of the vortex line(s) near 
the reconnecting point, which determines S{t), is still too small to satisfy this condition. 

The Biot-Savart model assumes scales much larger than the vortex core, which is effectively neglected. So it is 
not possible, as a matter of principle, to use the Biot-Savart model to study behaviour at the scale of the coherence 
length. Neither is possible, for practical computing reasons, to solve the GPE at the large scales explored by Biot- 
Savart calculations such as those in section (jlll F[) . 

The experimental observations of Paoletti et alS' agree with the Biot-Savart results in terms of the time symmetry 
and the exponent a = 0.5 of the power law (but it must be noticed that, for the GPE model, the average of the 
exponents before and after the reconnection is the same a sa 0.5 found in the experiment). In the experiment, the 
motion of the vortex lines was detected using solid hydrogen tracer particles of radius R 10~^ cm, ten thousand 
times larger than the coherence length in supcrfluid ''He. The rarefaction wave generated by GPE rcconncctions has 
a wavelength of about ten times the coherence length; although it is very deep initially, it quickly spreads out and 
vanishes as it moves away. On the scale of the tracer particles, the density of the fluid is thus constant, so it is not 
surprising that the Biot-Savart model is a better approximation to the observed dynamics of the vortex lines. 

Finally, the spread of the values of the coefficient A which we compute and the similar spread observed in the 
experiment is likely to arise from differences in the initial condition, geometry of nearby vortex lines and the velocity 
gradients which they induce, as discussed by Paoletti et alX. 
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Appendix A: Numerical method for the GPE 

Without loss of generality, and for sake of simplicity, we describe in detail the numerical method applied to the 
one-dimensional case of equation ^ and report at the end of the section the straightforward generalization to the 
three-dimensional case. It is convenient to split the GPE evolution in two time steps^: 

= ^VV, (Al) 

V^, = ^(1-|^P)^, (A2) 

thus separating linear and non-linear operators, where the subscript t denotes the time derivative. We assume that 
the solution is periodic in the domain a < x < i.e. i^{a) = ip{h), and seek a numerical solution in the time interval 
< i < T by expanding tj: via Fourier transform as 



'il){x,t)^ ^ (j)j{t)Fj{x), a<x<b, 0<t<T, (A3) 

■I 2 

where M is the number of modes, (t>j(t) arc the time-dependent Fourier coefficients, and the functions J^j{x) = 
-^^^q^'^^3t=^ are orthonormal (i.e. jj^TjTkdx = jj^ TjT-kdx = Sjk)- By computing the temporal and spatial 
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derivatives of t/i, substituting them in equation (|A2|) . multiplying the latter times J-^, integrating between a and b, 
and using the orthonormality property of the functions equation (|A2p reads 

^'k{t) = \)^k4>k{t), -Y<fc<Y' o<t<r, (A4) 



where \j = — is real and negative. The solution is, trivially, 



</,,(i) = e*i^''0fe(O). (A5) 
Since both t and A are real, = |(/'j(0)P for all j. Therefore, the total mass m is preserved: 

Ht)^\\m\\h = f H^.m^,t)dx^ E i0.-(Mi' = ii'^(o)iii. = m(o). (A6) 

J" 2 J~ 2 

Moreover, if (/) denotes the vector of Fourier coefficients (j) — [0„ m , ■ • • , then 

wmwh = mm = mom = m^) = m, m 

and the total mass m can be retrieved simply as the square of the norm of the complex-coefficient vector (p. 

The second part of the time-splitting scheme preserves mass as well. This is easy proved by taking the conjugate of 
equation \K2\ : from 4't = ~ we obtain ijj^ = (— i/2)(l — hence the derivative of with respect 

to t is 



2-(l - IV^I')^ 



0. (A8) 



dt dt ■'--r ■ -r-rt ^ 

The time independence of ji/'p is crucial because it implies that the solution of equation (jA2p is simply 

xP{x, t) = e*5(i-l'A(-.o)l')^(x, 0). (A9) 

The above method can be generalised naturally to the three-dimensional case. The unknown function ip(x,y, z,t) 
is expanded as 

^{x,y,z,t) = Y,Mt)^ji:x)J'k{y)J'i{z), (AlO) 

^ly 

where we use the notation Y., k i =Y. ■1_m^Y.,1_ My J2,1^m^- 

' J— 2 k — ^ t— 2 

After computing the temporal and spatial derivatives of ip, substituting them in equation (jA2p . multiplying the 
differential equation times Tm{x)J-'niy)J^s{z), integrating in space, and using the orthonormality property of the 
functions J-, the first part of the splitting now becomes 

(b'M*) ^ ^^Mjkiit), (All) 

where 



Jkl 



2TTj V f 27Tk V / 27rZ " ^ 



bx - J \by - ay J xb^- 



is real and negative and the mass-preserving solution is 
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As explained earlier, the second part of the splitting preserves mass too, and its solution is explicit. 

In order to outline the numerical algorithm it is convenient to formulate the partial differential equation ^ as an 
ordinary differential equation suppressing the spatial dependence and replacing '(/'('lO with u{t)^. We obtain the 
initial value problem 

u'{t) = [A + B{u{t))] u{t), u{to) = uo, (A12) 

Using second-order Strang splitting^, the solution u{t) — u(kAt) — Uk can be recursively determined by the multi- 
plication 



In our case A = ^^"^ ^(""(0) — 2^"*^ ^ previously shown, — \ip\'^ is constant, thus both 

operators A and B are linear. 

Since both solutions of the two parts of the splitting are explicit, the only numerical error introduced by the method 
is confined to the computation of the Fourier transform and its inverse. The second order error in time due to the 
Strang splitting can be improved to fourth order without further numerical complications^. 

In conclusion, assuming that the initial condition 'ip{x,y, z,to) is periodic in all spatial directions, the resulting 
algorithm is; 

1. y, z, io) = e^i'^~l'''''*°'^-'l ^ip{x,y, z,to): operator i? is applied in physical space 

2. ^(^o) = FFT('0(a;, J/, z, <o)): Fourier transform is applied 

3. (fijkiito + ^) = (^^^^^^''^<t'jki{to) V«, j, k: operator A = is applied in Fourier space 

4. il){x, y, z, to + ^) = IFFT((/)(io + ^))- inverse Fourier transform is applied to go back to physical space 

5. tp{x,y,z,to + At) = e^i^^-\'l'^''-y-''-^°+^^\^^{x,y,z,to + ^): operator B = -(l-|u(t)p) is applied in physical 



space 
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Clearly, step 5 is needed only to retrieve the physical solution ip^ix^y, z) at a certain time t ~ kAt, otherwise it can 
be avoided by merging step 5 and step 1 a single time step At. 

The limit of this scheme is that the initial condition ipoi^, y, z) must be periodic. If it is not, it must be made periodic 
by adding image vortices on a larger domain, i.e. more computational effort is required. However, certain geometries 
(such as two anti-parallel vortices aligned along the z-direction, moving along x and centered in {xol^yo)) allow us 
to double the grid points only in the z-direction due to the symmetry with respect to the y = and z = planes. 
In general, if all vortices are aligned along one axis (typically z), image vortices must be introduced in both other 
directions. For general geometries images vortices must be introduced in all directions, causing a memory allocation 
eight times larger. 
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angle /3 


^before 


Q^before 


-^after 


Oafter 


|^ = 7r 


1.36 


0.30 


1.88 


0.66 


l'^ = 


1.23 


0.39 


2.71 


0.63 




1.44 


0.41 


1.69 


0.68 


1^ = 1^ 


1.35 


0.44 


1.30 


0.69 


S'^ = 2 


1.41 


0.36 


1.01 


0.67 


h = h 


0.94 


0.42 


0.66 


0.73 


average 


1.29 


0.39 


1.54 


0.68 



TABLE I. Coefficients of the fit S{t) = A\t — to\" of the minimum distance between vortices before and after the reconnection 
at different initial angles /3. 
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FIG. 1. Left: Relative change of the dimensionless Hamiltonian energy E{t) as a function of dimensionless time t with 
respect to the initial energy Eo = E{0) during antiparallel reconnection. Right: Dimensionless density p as a function of the 
dimensionless radial coordinate r computed numerically (solid line, obtained by setting roo = 20 and using 500 grid points) 
and with the Fade approximation (empty circles). 




FIG. 2. Initial condition. Left: P = 37r/4. Right: /3 = -7r/2. 




FIG. 3. Snapshots of the evolution of two anti-parallel vortices (angle between vortices /3 
enhance the Crow instability, at i = (top left), t = 20 (top right), t = 30 (bottom left), t 
p = 0.2 are plotted to visualise the vortex cores. 



— tt), initially slightly perturbed to 

— 40 (bottom right). Isosurfaces of 
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FIG. 4. Minimum distance 5 between reconnecting vortices as a function of time t for the pair of reconnecting anti-parallel 
vortices (/? — tt) shown in figure Q. 




FIG. 5. Snapshots of the evolution of two initially straight vortices forming an angle /3 = 37r/4att = (top left), t — 
right), t = 20 (bottom left), t = 25 (bottom right). Isosurfaces of p = 0.2 are plotted to visualise the vortex cores. 
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FIG. 6. Snapshots of the evolution of two perpendicular vortices (angle between vortices /3 = 7r/2) ati = (top left), t 
(top right), t = 30 (bottom left), t — 40 (bottom right). Isosurfaces of p = 0.2 are plotted to visualise the vortex cores. 



= 20 
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t ~ to\ 1^ ~ ^0 



FIG. 7. Distance between vortices as a function of \t — to] before (top) and after (bottom) the reconnection for different values 
of the angle /3 between initial vortex lines (/3 = tt refers to anti-parallel vortices, and /3 — ■k/2 to orthogonal vortices). The 
computed values are joined by lines to guide the eye. The black solid lines are fits of the form S{t) = A\t — to|"; the fitting 
coefficients A and a are reported in table I. 
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FIG. 8. Isosurfaces ai p — 0.94 before (left, t = 37) and after (right and bottom, t — 54) reconnection for two anti-parallel 
vortices (angle between vortices (3 = tv). Note the mushroom-shaped rarefaction wave which moves away. 
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30.Q30.0 



FIG. 9. Isosurfaces at p and f = 54 as in figure ((Sj (bottom, anti-parallel vortices) to visualise the vortex core with superimposed 
the profile of p on the y = plane. Note the rarefaction wave. 
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FIG. 10. Snapshots of the evolution of the mushroom-shaped rarefaction wave ejected after reconnection, /3 = 7tt/8, at f = 53 
(top left), t — 55 (top right), t = 57 (bottom left), t = 59 (bottom right); isosurfaces of p = 0.94 are plotted to visualize the 
vortex cores and the wave as in figure (|8}. 
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FIG. 11. Snapshots of the evolution of two antiparallel vortices (angle l3 — tt) initially slightly perturbed to enhance the Crow 
instability, at t = 40 (top left), t = 80 (top right), t = 120 (bottom left) and t = 160 (bottom right). Isosurface of p = 0.2 are 
plotted to visualize the vortex cores. The initial condition is the same as in figure [3] but the computational box was extended 
to —60 < 2 < 60 to visualize the formation of vortex rings. 
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FIG. 12. Reconnection of two vortex rings initially set parallel to each other (as in the work of de Waele & Aarts^), computed 
with the Biot-Savart law: Top: at time t = 0; middle: to — t — 0.001 s; bottom: t — to = 0.005 s. The vortex lines are 
colour-coded to indicate the magnitude of the velocity (in cm/s, see legend on each figure). The box is for visualization only. 
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FIG. 13. Reconnection of two vortex rings initially set perpendicular to each other, computed with the Biot-Savart law: Top: 
at time t — 0; middle: to — t — 0.001 s; bottom: t — to — 0.005 s. The vortex lines are colour-coded to indicate the magnitude 
of the velocity (in cm/s, see legend on each figure). The box is for visualization only. 
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FIG. 14. Minimum distance between the filaments 5{t) (cm) vs {to — t) (s) before the reconnection (top) and vs {t — to) 
after the reconnection (bottom), corresponding to the Biot-Savart evolution of two parallel vortex rings shown in figure (lllf) 
(solid black line) and of two perpendicular vortex rings shown in figure (|12|) (dot-dashed blue line), to is the time at which 
the reconnection takes place. The dashed red line expresses 5{t) = nito — t)/(27r) found by de Waele & Aarts^i (top) and 
&{t) = ^J-KK[t — to) (bottom) as a guide to the eye. 



